Energy loss for heavy quarks in relation to light partons; is radiative energy loss for 

heavy quarks anomalous? 
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The scaling properties of jet suppression measurements are compared for non-photonic electrons 
(e ) and neutral pions (ir°) in Au + Au collisions at y^ijvjv = 200 GeV. For a broad range of trans- 
verse momenta and collision centralities, the comparison is consistent with jet quenching dominated 
by radiative energy loss for both heavy and light partons. Less quenching is indicated for heavy 
quarks via e ; this gives an independent estimate of the transport coefficient q that agrees with its 
magnitude obtained from quenching of light partons via 7r°'s. 



Partonic energy loss provides an important tomo- 
graphic probe of the hot and dense matter produced 
in high energy nuclear collisions Such energy loss 
is manifested as a suppression of high-px hadron yields 
(jet quenching) in central and midccntral A+A collisions, 
when compared to the binary-scaled yields from p+p col- 
lisions 0, Q . There is mounting evidence that the pri- 
mary mechanism for light-parton energy loss is medium- 
induced gluon bremsstrahlung in the quark-gluon plasma 
(QGP) produced in A+A collision Q, [1, 1, @, 0, i, H E3] • 
Therefore, measurements of jet quenching can provide 
important constraints for determining the transport co- 
efficients of the QGP. 

Further insight into properties of the QGP can be 
gained from the production and propagation of particles 
carrying heavy quarks (charm or bottom). However, a 
central issue germane to the development of this probe, 
is the unsettled question of the quantitative difference be- 
tween the quenching for light and heavy partons in hot 
QCD matter. Dokshitzcr and Kharzeev |7fl have argued 
that much less quenching is to be expected from heavy 
quarks because their associated gluon radiation is sup- 
pressed by the dead cone effect. In contrast, Zakharov 
[ill . fl2| argues that the radiative energy loss for a finite- 
size QGP could even have an anomalous mass depen- 
dence (ie. energy loss for heavy quarks greater than that 
for light quarks) which stem from quantum final state 
effects when the gluon formation time is comparable to 
the "size" of the QGP. 

Here, we use the published jet quenching measure- 
ments for non-photonic electrons (e^) (l3j and neutral 
pions (ir°) [3, [l5[ at y/s^N = 200 GeV to compare, as 
well as to quantify the difference in the quenching pat- 
terns of the transverse momentum spectra for jets pro- 
duced from scattered light- and heavy partons. We then 
use this difference to constrain an estimate of the trans- 
port coefficient q and the ratio of viscosity to entropy 
density rj/s. We find that the value of these new es- 
timates (of q and rj/s) are similar to the ones obtained 
via jet quenching measurements for light partons via 7r°'s 

The nuclear modification factor (Raa) is used to quan- 



tify the magnitude of jet suppression in A+A collisions 

1; 
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where a pp is the particle production cross section in p+p 
collisions and (Taa) is the nuclear thickness function av- 
eraged over the impact parameter range associated with 
a given centrality selection 
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The corresponding average number of nucleon-nucleon 
collisions, (N co u) — cr™ 6 ' (Taa), is obtained with a 
Monte-Carlo Glauber-based model calculation (Itl Il7j]. 

To facilitate our comparisons of the quenching patterns 
of the transverse momentum spectra for jets produced 
from scattered light (/)- and heavy (h) partons, we exploit 
the formalism of Dokshitzer and Kharzeev (DK) Q; 
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where a s is the strong interaction coupling strength, C F 
is the color factor, L is the path length [of the medium] 
that the parton traverses, M is the mass of the heavy 
parton and q is the transport coefficient of the medium. 

The first term in the exponent in Eq. [T] represents the 
quenching of the transverse momentum spectrum which 
is the same for both light and heavy partons Q , 
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(modulo the difference of the £ parameters determined 
by the pt distributions in the vacuum) . The second term 
in the exponent is specific for heavy quarks and is a direct 
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consequence of the fact that the mass of the heavy quark 
leads to a suppression of gluon radiation. Hence the pre- 
diction that the suppression for the heavy hadron px dis- 
tribution should be smaller than that for light hadrons. 

The magnitude of the heavy-to- light suppression factor 
can also be estimated as Q; 



Raa(pt,L) 
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which clearly depends on the size (L) of the QGP medium 
and the properties encoded in the value of its transport 
coefficient q. 

An essential point reflected in Eqs. Q] - [3] is that they 
give specific testable scaling predictions for the px and L 
dependence of R AA (px,L), and the heavy-to- light sup- 

' should scale 



r aa(pt,L) 



pression ratio [h}. That is, In 

as L and l/y/px respectively, and the heavy-to- light sup- 
pression ratio should scale as L. As discussed below, 
these scaling properties provide crucial validation tests 
of our analysis framework, as well as an important con- 
straint for estimating q and rj/s. 

The result of one such test for Cu+Cu and Au+Au 
collisions is given in Fig. [TJ where we show a plot of 
In [Raa(pt, L)\ vs. 1/ v / Pt for 7r°'s. Here, results are 
shown for px > 5 GeV/c in panel (a) and for 2.75 < 
Pt ^5 5 GcV/c in panel (b). Both panels indicate the 
predicted linear dependence on l/^/px (cf. Eq.[2]), albeit 
with different slopes for each px range. A comparison 
of the results for the two collision systems also indicate 
similar quenching magnitudes when the size of the re- 
spective collision medium is comparable. Note that bet- 
ter agreement is achieved if the small difference in the 
length estimate arising from the respective centrality cut 
is taken into account. 

The inset in Fig.[IJb) confirms the predicted linear de- 
pendence of In [i?AA (.Pt i L)] on L for Cu+Cu collisions. A 
similar de pen dence has been reported earlier for Au+Au 
collisions [id ]. Here, as in Ref. we have used the 

transverse size of the system R as an estimate of the oper- 
ative path length L. To obtain this estimate, the number 
of participant nucleons iV part was estimated for each cen- 
trality selection, via a Monte-Carlo Glauber-based model 
fl6l , 03] ■ The corresponding transverse size R was then 
determined from the distribution of these nucleons in the 
transverse {x, y) plane via the same Monte-Carlo Glauber 



model: ^/R = y y-^r + -^r) > where a x and a y are the re- 
spective root-mean-square widths of the density distribu- 
tions. For these evaluations, the initial entropy profile in 
the transverse plane was assumed to be proportional to a 
linear combination of the number density of participants 
and binary collisions [H, O • The latter assures that the 
entropy density weighting is constrained by multiplicity 
measurements. 

Figure [5] shows a comparison of the Raa measurements 
obtained for neutral pions (w ) and non-photonic elec- 
trons (e^) in minimum-bias Au+Au collisions (T3 , [Til ]. 
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FIG. 1: ln[R AA {pT,L)] vs. l/-y/p¥ for p T > 5 GeV/c (a) 
and 2.75 < p T < 5 GeV/c (b) for 7r°'s produced in Cu+Cu 
(N part = 99) and Au+Au (N part = 114) collisions. The inset 
shows In [Raa(pt, L)\ vs. L for tv°'s with pri = 9.25 GeV/c 
and pT2 = 5.75 GeV/c in Cu+Cu collisions. Error bars are 
statistical only. The curves in each panel are linear fits to 
each data set (see text). 



As in Fig. [T] results arc shown for px > 5 GeV/c in 
panel (a) and for 2.75 < px < 5 GeV/c in panel (b). 
The latter px range for is dominated by decay con- 
tributions from D-mesons. However, as demonstrated in 
recent measurements pol [2~i| , there is an increasing role 
of the B-meson contributions to the non-photonic elecron 
spectrum for px > 5 GeV/c. Such contributions could 
serve to complicate the interpretation of the Raa mea- 
surements for above px ~ 5 GcV/c. 

The 7T° measurements in Fig. Ufa) show the predicted 
linear dependence on 1/^/px. However, the statistical 
significance of the data for e^- does not allow a similarly 
decisive conclusion in this px range. The measurements 
shown in Fig.[2^b) contrast with those of Fig. [2ja) . Here, 
both data sets confirm the the predicted linear depen- 
dence on 1/yfpx. Even more striking is the observation 
that the magnitude of the quenching for 7r° is signifi- 
cantly larger than that for e ± . This observation is in 
accord with the prediction of Dokshitzer and Kharzeev 
(cf. Eq. Q] and Ref. 0) that the suppression of heavy 
hadrons should be less than that for light hadrons. As 
discussed below, it is also congruent with the observa- 
tion that the magnitude of the azimuthal anisotropy for 
e (characterized by the second Fourier coefficient v?) is 
much less than that for pions (for p T > 2.5 GeV/c) jl3| j . 

Figure [3] compares and contrasts the L dependence of 
In [Raa ijpT for tt° and e ± for the px selections indi- 
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FIG. 2: Iii[Raa(Pt,L)] vs. for p T > 5 GeV/c (a) 

and 2.75 < pr < 5 GeV/c (b) for minimum-bias Au+Au 
collisions. Error bars are statistical only. The curves in each 
panel are linear fits to each data set (see text). 
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FIG. 3: (a) In [Raa(pt, L)] vs. L for n and non-photonic 
electrons (e ) as indicated. The dashed and dot-dashed 
curves indicates a linear fit to each data set. The error bars 
in both panels are statistical only, (b) In [( ^f" A | PT '^j vs. 

L. This ratio is obtained from the fits to the it and e data 
sets in panel (a). The dotted lines represent the statistical 
error band. 



cated. We emphasize here that kinematic studies show 
that the pr selection for e ± (cf. Fig. [3|) leads to a domi- 
nant decay contribution from D-mesons of higher px- 
The dashed and dot-dashed curves in panel (a) are lin- 
ear fits to the data set for it and respectively (the 
data point for the most peripheral collision is treated as 
an outlier). The dashed curve in panel (b) shows the L 



dependence of the ratio In 



obtained from 



?aa(pt,L) ' 
Aa(Pt,L) / 

the fits in panel (a). These curves validate the predicted 
linear dependence on L (cf. Eqs. [T]-[2]), as well as the 
attendant difference in slope for tt° and e^. We interpret 
the latter as an independent indication that the mecha- 
nism for jet quenching is dominated by radiative energy 
loss. 

The fits to the data in Fig. [3] (a) indicate intercepts 
of L = 0.6 ± 0.1 fm and L = 0.9 ± 0.15 fm for the light- 
and heavy mesons respectively. This suggests a minimum 
path length requirement for the initiation of jet suppres- 
sion for both light and heavy partons, but with the pos- 
sibility of a larger path length requirement for the heavy 



quark. The corresponding ratio In ( ^ AA jp T '^j J ; plot- 
ted as a function of L in Fig.[3](b), yields a slope value of 
0.4 ± 0.04 fm -1 . This value reflects the difference in slope 
for 7T° and shown in Fig. [3ja). It is noteworthy that 
the latter is fully compatible with the observed difference 
of about a factor of two in the measured magnitude of 
v 2 for 7T° and e ± for p T > 2.5 GeV/c [HI- Thus, it pro- 
vides further confirmation that the azimuthal anisotropy 
of particle yields (for the meson px ranges of interest) 



stem from the path-length dependence of jet quenching. 
This observation further underscores the importance of 
reporting jet suppression measurements in conjunction 
with anisotropy (i^) measurements at high px- 

To obtain an "independent" estimate of the magni- 
tude of q, we assume a negligible e ± contribution from 
B-mesons and use the slope (0.4 ± 0.04 fm -1 ) extracted 
from Fig. 0(b) in concert with Eq. [31 This gives the 
value q w 0.73 ± 0.12 GeV 2 /fm for the values a s = 0.3 
[2|, C F = 4/3, M = 1.5 GcV and (p T ) ss 6 GeV/c (for 
D-mesons). This estimate of q is comparable to the re- 
cent estimate of ~ 0.75 GeV 2 /fm obtained solely from 
7T° suppression measurements with the same definition 
of L [H|. It is also compatible with the value of 1 - 2 
GeV 2 /fm obtained from fits to hadron suppression data 
within the framework of the twist expansion [23|, [24[ . 

The value for q can be used to estimate the ratio of 
the shear viscosity (77) to the entropy density (s) as (25|; 



1.25^ 



T 



where T is the temperature. This gives the estimate 
47rf w 1.4 ± 0.4 for T - 220 MeV [H, which is com- 
parable to the value estimated via ir° suppression mea- 
surements [Io| and flow measurements [ljl ■ We con- 
clude therefore, that the relatively short mean free path 
in the QGP [lj| leads to hydrodynamic-like flow with 
small shear viscosity, as well as significant quenching of 
both light and heavy partons. 
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In summary, we have compared the scaling proper- 
tics of jet quenching for both light- and heavy partons, 
via suppression measurements for tt° and e^. Our com- 
parisons confirm the predicted pr and L dependence for 
medium-induced gluon radiation in a hot and dense QGP. 
The difference in the magnitude of the quenching for 
light- and heavy partons, do not give a strong indication 
for anomalous heavy quark energy loss in the pr range 
studied. Instead, it is compatible with the prediction 
that gluon radiation from heavy quarks is suppressed by 
the dead cone effect. This difference also gives an esti- 
mate q ~ 1 GeV 2 /fm, which is essentially the same as 
that obtained via ir° suppression measurements. Future 
detailed measurements of D- and B- mesons at high pr, 
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as well as model calculations which take full account of 
the reaction dynamics, will undoubtedly provide more 
detailed mechanistic insights for heavy quark energy loss 
and improved constraints for the transport properties of 
the QGP. 
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